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Abstract 

We address a recently introduced model describing a system of periodically 
coupled nonlinear phase oscillators submitted to multiplicative white noises, 
wherein a ratchet-like transport mechanism arises through a symmetry- 
breaking noise-induced nonequilibrium phase transition. Numerical simu- 
lations of this system reveal amazing novel features such as negative zero- 
bias conductance and anomalous hysteresis, explained resorting to a strong- 
coupling analysis in the thermodynamic limit. Using an explicit mean-field 
approximation we explore the whole ordered phase finding a transition from 
anomalous to normal hysteresis inside this phase, estimating its locus and 
identifying (within this scheme) a mechanism whereby it takes place. 
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I. INTRODUCTION 



Feynman's ratchet-and-pawl example , illustrating the impossibility for a microscopic 
rectifying device to extract work in a cyclic manner from the equilibrium fluctuations of 
a single heat bath, spurred in turn the search for heat engines operating in a far from 
equilibrium regime between two heat baths. The field of "nanomechanics" (more specifically, 
that of noise-induced transport or "Brownian motors") is now about one decade old M. In 
the early works, a requisite for these devices to operate (besides their obvious built-in, 
ratchet-like, bias) seemed to be that the fluctuations be correlated ||. That requirement 
was relaxed when "pulsating" ratchets were discovered: in these it is the random switching 
between uncorrelated noise sources which is responsible of the rectifying effect [|J. 

A recent new twist has been to relax also the requirement of a built-in bias || : a system 
of periodically coupled nonlinear phase oscillators in a symmetric "pulsating" environment 
has been shown to undergo a noise-induced nonequilibrium phase transition, wherein the 
spontaneous symmetry breakdown of the stationary probability distribution gives rise to 
an effective ratchet-like potential. The authors introduced the aforementioned mechanism 
and its striking consequences, such as the appearance of negative zero-bias conductance and 
anomalous hysteresis, which they illustrated through numerical simulations and explained 
by resort to the strong-coupling limit. By anomalous hysteresis we refer to the case where 
the cycle runs clockwise, in opposition to the normal one (as typified by a ferromagnet) that 
runs counterclockwise. 

Exploiting our previous experience in a lattice model displaying (like the present one) 
a symmetry-breaking nonequilibrium phase transition || and in order to set a firm ground 
for further work, we addressed the model using an explicit mean-field approach [0], focusing 
on the relationship between the shape of the stationary probability distribution (as well as 
the number of solutions to the mean-field equations) and the transport properties in its 
different regions. Hence it is our aim in this work to report on the thorough exploration of 
the ordered phase, on the characterization of its subregions, and on features related to the 
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transition from anomalous to normal hysteresis in the behavior of the particle current as a 
function of the bias force. Our main finding is that there exists a close relationship between 
the character of the hysteresis loop on one hand, and the shape of the stationary probability 
distribution as well as the number of "homogeneous" solutions (a term to be clarified later) 
on the other. 

In the following sections we successively introduce the model, describe the mean-field 
approach, discuss our numerical results and draw our conclusions. For the benefit of the 
reader we have included an appendix where some (in our judgment) subtle calculations are 
performed in some detail. 

II. THE MODEL 

In Ref. |J the authors consider a set of globally coupled stochastic equations of motion 
(in the overdamped regime, and to be interpreted in the sense of Stratonovich) for the N 
degrees of freedom (phases) Xi(t): 

Xi = -|| + V2fUt) - ^ E K(Xi - Xj). (1) 

The model just set up can be visualized (at least for some parameter values) as a set of 
overdamped interacting pendula. The second term models, as usual, the effect of thermal 
fluctuations: T is the temperature of the environment and the are additive Gaussian 
white noises with zero mean and variance one 

&(*)> = o, <&(*)&(0> = M(*-0- (2) 

The "pulsating" potentials Ui(x,t) are among the key ingredients in the model || (x G 
[-L/2, L/2] is a phaselike real variable that runs over the range of Xi(t), namely the allowed 
values for any realization of any X,i at any time t). They consist of a static part V(x) 
and a fluctuating one: Gaussian white noises f]i(t) with zero mean and variance one (i.e. 
obeying also Eq. (|2|) but assumed nonthermal in origin) are coupled multiplicatively (with 
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intensity Q) through a function W(x). Even though we adopted L = 2tt in all our numerical 
calculations, we kept in all the equations the dependence on L in order to have the most 
general expressions. For the analysis of the noise-induced ratchet effect a "load force" F, 
producing an additional bias, is also included 

Ui(x,t) = V(x) + W(x)^2Qni(t)-Fx. (3) 

As already stated, both V(x) and W(x) are assumed to be periodic (period L) and moreover 
they are symmetric: 

V(-x) = V(x) , W(-x) = W(x) 

which means that there is no built-in ratchet effect. In Ref. || the choice was 

V(x) = W(x) = — cos x — A cos 2x, (4) 

with A > so as to remove an accidental degeneracy hindering the spontaneous symmetry 
breakdown ||, but not strong enough to create a local minimum at L/2 = ir. With the choice 
A > the direction of the particle current (X) turns out to be opposite to that of symmetry 
breaking in the stationary probability distribution P st (x), and it is this effect which leads 
in turn to such oddities as negative zero-bias conductance and anomalous hysteresis || (a 
nice animation illustrating this phenomenon can be found on the web |§). 

The interaction force K(x — y) = —K(y — x) between oscillators (the other key ingredient 
in the model) is also assumed to be a periodic function of x — y with the same period L = 2tt 
as V(x), W(x). In Ref. § it is 

K(x)=K smx, K >0. (5) 

As indicated before, the model can be visualized (at least for A — » 0) as a set of overdamped 
and interacting pendula (only their phases matter, not their locations) interacting with one 
another through a force proportional to the sine of their phase difference (a force that is 
always attractive in the reduced interval — tt < x — y < n). 
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Summarizing, the complete set of parameters in the model is: N, F, T, A, Kq and Q. 
Except for numerical simulations the exact value of N is unimportant, as long as it is large. 
As already said, F is just an auxiliary tool for the analysis. We shall fix the values of T and 
A as in Ref. ||, namely T = 2,A = 0.15. So the important parameters in the model are 
Kq (governing mostly the "drift" terms in this set of generalized Langevin equations) and 
Q (governing mostly the "diffusion" ones). As discussed in Ref. ||], the Gaussian character 
of rji(t) allows it to be added to Hence it suffices to consider the ?7i(i)'s, now coupled 

through S(x) = V2[T + Q(W') 2 ] 1/2 (note that S(x) = V2g(x) as defined in Ref. (1). 

III. MEAN-FIELD ANALYSIS 

On account of the choice made in Eq. (£|), the interparticle interaction term in Eq. (|]) 
can be cast in the form: 

1 N 

- £ K(Xi - X,) = K [Ci(t) sinX, - Si(t) cosX t ] . (6) 

7V 3=1 

For N — ► oo we may approximate Eq. (|l|) a la Curie- Weiss, replacing Ci(t) = 
N' 1 J2j cos Xj(t) and Si(t) = N^ 1 J2j siia Xj(t) by C m = (cosxj) and S m = (siriXj) respec- 
tively, to be determined as usual by self-consistency. This decouples the system of stochastic 
differential equations (SDE) in Eq. ([I]), which reduces to essentially one Markovian SDE for 
the single stochastic process X(t): 

X = R(X) + S(X) V (t), (7) 

with 

R{x) = -V'(x) +F- K m {x) 

= — sinx(l + K C m + 4v4cosa;) + A'o'S'mCosx + F (8) 

(where K m (x) = K [C m sin x — S m cosx]) and 

S(x) = \/2{T + Q[W'(x)} 2 } = ^]2{T + g[sin x + 2 A sin 2a;] 2 }, (9) 
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so that the "spurious" contribution to the drift in the Stratonovich interpretation is 

-S(x)S'(x) = QW'{x)W"(x) =g(sinx + 2Asin2x)(cosx + 4Acos2x). (10) 

We may write R(x) = —V'(x) in terms of an effective (since C m and S m are determined by 
self-consistency in terms of P st (x) below and moreover by the asymmetric, hence ratchet- 
like) potential 

V{x) = V{x)+ f X K m {y)dy-Fx 
Jo 

= — [cosx(l + K C m ) + Acos2x] — K S m smx — Fx. 
Similarly S(x)S'(x)/2 can be derived from — Q[W'(x)} 2 /2. 

A. The stationary probability distribution function 

The Fokker-Planck equation associated with the SDE in Eq. (0) is 

d t P{x,t) = d x {-[R{x) + ^S(x) S'{x)]P(x,t)} + l -d xx [S\x)P(x,t)} (11) 

(see the appendix, Eqs. ( p5| ) and (^)) and its normalized stationary solution with periodic 
boundary conditions and current density J ^ is [fj.^[ 

p -4>(x) jj( t \ 



where 



<*> = - 2 f*lsr- (13) 



H ^=f L ^-^r- < 14 > 

and JV = j_L/ 2 dxP st (x). The positive sign of S(x) and the exponentials implies that of 
H(x) and hence that of P st (x) and J\f, as it should be. On the other hand, although R(x), 
S(x) and P st (x) are periodic by construction, <p(x) is not required to be so; in fact it increases 
on each cycle by an amount 
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1 r L dxp(x) sin x K S m r L dx cos x F r L dx , . 

0( > = fJo q^) T~ Jo ~^T^TJo WY ' ^ 

with p{x) = 1 + KoC m + 4Acosx and q(x) = 1 + (Q/T) sin 2 a; (1 + 4Acosx) 2 , both even 
functions of x. Since only the first term vanishes identically, for nonzero F or S m it will 
be the case generically that (f>(L) ^ 0; thus the form of H(x) in Eq. (|T3]) is designed to 
compensate for this fact. Moreover (as shown in the appendix) for A > 0, [g(x)] _1 (> 0) gives 
less weight to the positive cos x values than to the negative ones, and it does increasingly so 
the larger Q is; hence the two nonzero contributions to 4>(L) compete with each other. 
According to Eq. (|37|) in the appendix 

J = [1 - e ^ (L) ]/2A/-, (16) 

hence the sign of J is that of 1 — e^ L ^ and — on the other hand — the "holonomy" condition 
e <t>{L) _ j j m plies J = and H(x) = const, (see the appendix, Eq. (|30D). Equation ( |16D is 
a self-consistency relation since both J\f and 4>(L) keep information on the shape of P st (x) 
(in the latter case through C m and S m ). A nonzero J is always associated with a symmetry 
breakdown in P st (x) (namely, P st (—x) ^ P st (x)). This may be either spontaneous (our 
main concern here) or induced by a nonzero F. 

B. The self-consistency equations 

As indicated earlier, the stationary probability distribution P st (x) depends on both S m 
and C m , since R(x) includes K m (x). Their values arise from requiring self-consistency, which 
amounts to solving the following system of nonlinear integral equations 

pL/2 



F cm = C m , with F cm = (cosx) = dx cosx P st (x, C m , S m ), (17) 

J -L/2 
[L/2 

F sm = S m , with F sm = (sin a;) = / dx sinx P st (x,C m , S m ). (18) 

J -L/2 

These equations give C m and S m for each set of the parameters (Q, K ) that define the state 
of the system (assuming T, A and F fixed). 



7 



For F = 0, the choice S m = makes R(x) in Eq. (Q) an odd function of x; this in turn 
makes <fi{x) in Eq. ( |i3|) even, and then the periodicity of P st (x) in Eq. ( |T2| ) (in the form 
P st (—x) = P st (—x — L)) implies that the stationary probability distribution is also an even 
function of x. So the problem of self-consistency reduces to the numerical search of solutions 
of Eq. (PTT|) , with S m = 0. Nonetheless, a plausibility argument leads to an intuition on the 
existence of some solutions of this integral equation (and their stability) in this symmetric 
case: 

• Since — as argued — J = and H(x) = const, holds, it turns out from Eq. (|Tl|) that the 
set of critical points x c of P st (x) must obey R(x c ) = \S(x c ) S'(x c ). 

• Since x c = 0,ir (mod 2ir) belong to that set, a possible way to satisfy the integral 
equation in Eq. ( |17D is for P st (x) be concentrated around those values. 

• Then, by analogy with a pendulum, one expects the solution with C m > to be the 
stable one for K Q /Q large enough (what we shall call an "interaction-driven regime" 
or idr). For K /Q small enough ("noise-driven regime" or ndr) the stable solution 
can have C m < (implying an "angle" larger than 7r/2 with respect to the mean field) 
since it corresponds to shaking the pendula violently. 

Beyond these handwaving arguments it must be said that, since cosx in Eq. (|T7|) is an even 
function of x, in order to determine the stability of the true solutions it suffices to use the 
Curie- Weiss (i.e., the one-parameter) criterion, namely to check whether the slope at S m of 
the integral in Eq. (|18D is less or greater than one. As a complementary check, a small-x 
expansion of <f)(x) (see the appendix) confirms that P st (x) is indeed Gaussian at x — 0. For 
small F ^ 0, P st (x) gets multiplied (in this approximation) by exp [Fx/T\{?z 1 + Fx/T) 
which leads to a nonzero value of S m = kF, with k > 0. By the mechanism discussed after 
Eq. (|To]), for Q large enough <j)(L) > and by Eq. flToT) J < 0. As will be shown later, this 
effect manifests itself in a negative zero-bias conductance. 

We conclude that, for F = there are always one or more solutions to Eqs. (fTT|) — (P~5|) 
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with S m = and one of these is the stable one in the "disordered" phase. As argued in Ref. 
0, for N — > oo a noise- induced nonequilibrium transition takes place generically towards 
an "ordered" phase where P st (—x) 7^ P st (x). In the present scheme this asymmetry should 
be evidenced by the fact that the solution with S m = becomes unstable in favor of other 
two solutions such that -P 2 st ( x ) = Pfi. - x )i characterized by nonzero values ±|£ m |. This fact 
confers S m the character of an order parameter. 



C. The phase boundary 



Since sin a; is an antisymmetric function, Eq. flI5|) results impractical for the task of 
finding the curve that separates the ordered phase from the disordered one, given that on 
that curve S m is still zero. For that goal (exclusively) we solve, instead of Eqs. flT7|)- 
the following system: 

-L/2 

L/2 dpst 

ax sinx — — = 1. (20) 



/•L/2 

J l ^dx cosxP st {x,C m ,0) = C m , (19) 



L/2 3, 



S' 



Sm=0 



D. The particle current 

The appearance of a ratchet effect amounts to the existence of a nonvanishing drift 
term (X) in the stationary state, in the absence of any forcing (F = 0); in other words, 
the pendula become rotators in an average sense. As discussed above, the cause of this 
spontaneous particle current is the noise-induced asymmetry in P st (x) ||. 

As it is shown in the appendix 

(X) = f L/2 dx \r(x) + \s{x)S'{x)} P st (x, C m , S m ), (21) 

J -L/2 L 2 



and the final result is 



f 1 _ e <KL) 1 

»="=hr- \ L - (22) 



Hence (X) has the sign of J and can be also regarded as an order parameter. In fact, 
from Eqs. (§) and ( |21|) — or equivalently from Eqs. (§), (0) an d fl22|) — one may suspect the 
existence of a tight relationship between (X) and S m . 

IV. NUMERICAL RESULTS 

Figure [l] displays (on the same scale as Fig. lb of Ref . || , with which it fully coincides) 
the phase diagram obtained by solving Eqs. (|l9l) -(f20"l) using the Newton-Raphson method. 
In the region above the full line ( "ordered region" ) the stable solution to Eqs. (|l7|)- ([r^) has 
S m 7^ 0. Notice that this noise-induced phase transition is reentrant: as Q increases for 
K Q = const., the "disordered phase" (S m = 0) is met again. The multiplicity of mean-field 
solutions in the ordered region, together with the fact that some of them may suddenly 
disappear as either K or Q are varied (a fact that, as we shall see, is closely related to the 
occurrence of anomalous hysteresis) hinder the pick of the right solution in this region. 

A more systematic characterization of the aforementioned multiple solutions is achieved 
when the branch to which they belong is followed from its corresponding "homogeneous" 
(S m — 0) solution. Accordingly, the dashed line in Fig. [1| separates two sectors within the 
ordered region with regard to the homogeneous solutions. Below it ( "noise-driven regime" or 
ndr) there is a single solution with S m = and C m < (as already suggested, in this regime 
a solution with C m < can be stable since it corresponds to shaking violently the pendula). 
Above it ("interaction-driven regime" or idr) there are three: two of them have opposite 
signs and (for K Q /Q large enough) |C m | ~ 0.9; the remaining one has C m ~ 0. Note that 
this line presents a dip, whose meaning will be discussed in relation with the character of the 
hysteresis loop. We have studied the shape of P st (x) and the behavior of (X) as a function 
of F for different locations in this (Q, K ) diagram. The squares in Fig. [I] indicate several 
positions inside and outside the ordered zone for K = 10. For this value, the separatrix 
between both regimes lies around Q = 6. Figure |2| illustrates the crossing of the dashed line 
in Fig. 1, and the persistence of the negative solution after the disappearance of the other 
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two. 



A. Analysis at constant coupling 

Figure |3| shows (for the true solution, namely the stable S m 7^ one) the evolution of 
P clS <X function of Q, for Kq = 10: whereas for Q = 1 and Q = 21 it is a symmetric 
function of x, there is a spontaneous breakdown of this symmetry for the remaining values 
(the system has to choose between two possible asymmetric solutions, of which just one is 
shown). A noticeable feature of P st (x) is that for some Q > 6 it becomes bimodal (in fact, 
already for Q = 6 do we see an indication that a second peak is developing). According to 
Fig. H| and Eq. (|17|) , for Q < 6 it may be expected that (for the true solution) C m be positive 
and even relatively large (however, between Q = 2 and Q = 6 other solutions are possible, 
which lead to other shapes of P st (x) not shown). As the second peak develops and becomes 
higher than the original one, C m shifts toward small negative values (in this region, the one 
depicted is the only possible solution). 

The squares in Fig. |4] plot (always for Kq = 10) the two possible values of the spontaneous 
drift velocity (X) ^ in the ordered phase as functions of Q. The vertical thick line indicates 
the value of Q at which the transition from anomalous to normal hysteresis occurs for 
Kq = 10. The effect of a moderate positive bias F on (X) in the normal region (at the 
right of the thick line) is clearly understandable (the only surprising feature is that on the 
reentrant branch of the phase boundary the transition for F 7^ is so steep that it resembles 
a first order one). But the most striking feature is the sudden disappearance of a "forward" 
particle current for the shown values of F as we cross the thick line towards the left (although 
it may still exist for lower values of F): this is the manifestation of the anomalous hysteresis. 
As suggested in the first paragraph of this section, this phenomenon is intimately related to 
the sudden disappearance of some of the multiple solutions when either Kq or Q are varied 
in the idr (see Fig. 0). 

Figures |5|(a) to |5|(d) present a sequence of (X) vs F plots, varying Q across the thick line 
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of Fig. £2[ In these, all the solutions to Eqs. (|T7| ) -([T8| ) but the one belonging to the branch 
starting at C m ~ for S m = have been included. For Q = 5.97 (Fig. ||(a)) two (unstable) 
solutions meet at (X) = for F — 0. The progressive withdrawal of one of them out of 
the F ~ region with increasing Q until its complete disappearance (figs. ^|(b) to |](d)) can 
be traced back (through their corresponding branches) to the disappearance of solutions for 
S m = 0. Moreover, it is only after this solution has completely disappeared that the stable 
solution begins to exist for larger values of F and thus normal hysteresis sets in (Fig. ||](d)). 

It is also instructive to see how the different branches in figs. ||(a) to f|(c) develop as one 
enters the ordered region from the left. Figure |6] shows for three points on the Kq = 10 
line the (X) vs F plots for all the existing solutions (but not the one belonging to the 
branch starting at C m ~ for S m = 0). For Q = 1 there is a single stable solution 
displaying negative zero-bias conductance; for Q = 1.7 (right on the phase boundary) a 
second (unstable) solution appears and the anomalous hysteretic behavior (clearly seen for 
Q — 3) sets in. As suggested by Fig. f|, the situation is different at the reentry: Figure 
shows that the disappearance of the (normal) hysteretic behavior at the phase boundary in 
the ndr is in fact abrupt, signalling a first order phase transition in this regime. 

B. Analysis at constant noise intensity 

A complementary view of the transition in figs. ||(a) to §(d) is obtained by varying Kq 
at Q = 6.0 (figs. ||(a) and |8](b)): for Kq = 7.25 a very small normal hysteresis loop can 
be appreciated, which has grown rather large already for Kq = 8.0; for Kq = 10.0 the loop 
has become anomalous and a third branch has appeared, forming a cusp at the endpoints 
of the loop; for a larger Kq (Fig. ||(b)) the cusp develops into a curl. For larger values of 
Q (always across the dashed line in Fig. [I]) the general pattern is about the same (see Fig. 
|9] for Q = 9.5, where the kink of the normal loop at a position rather close to the endpoint 
of the anomalous one is suggestive); a similar plot to that in Fig. |2] (but now varying Kq at 
Q = 10.0) is shown in Fig. [Tl], where the remaining solution displays a larger value of \C m \. 
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Finally, it is interesting to elucidate the nature of the transition at the left of the dip in the 
dashed line in Fig. [I]: as Fig. [11] shows, here the loss of two solutions is not accompanied by 
a change in the character of the hysteresis loop, which remains anomalous. 



V. CONCLUSIONS 

We have shown the existence of a sharp transition in the behavior of the system inside 
the ordered phase, from an "interaction-driven regime" (idr) (typically for K /Q larger than 
about three halves) towards a "noise-driven regime" (ndr) which differs from the former in 
several aspects: 

a. Although (X) shows hysteretic behavior as a function of F everywhere inside the 
ordered phase, in the idr its character is anomalous (namely, clockwise) whereas in the 
ndr it is normal (counterclockwise). Moreover, whereas the height of the anomalous 
hysteresis loop increases continuously at the phase boundary in the idr ((X) acts as an 
order parameter in a second order phase transition), the disappearance of the normal 
one proceeds by shrinking its width at a more or less finite height (the transition at 
the reentry is of second order but it is so steep that resembles a first order one). 

b. The shape of the stationary probability distribution function (pdf) changes qualita- 
tively in going from the idr to the ndr (it becomes bimodal and remains so as the 
disordered region is reentered and the pdf becomes symmetric again, the peak at it 
then being higher than the one at 0). 

c. Whereas in the idr there are several solutions with S m = to the mean-field equations 
(Eqs. ([17| ) -([T8| ) ), in the ndr there is a unique solution with S m = 0. Solutions of this 



kind are relevant as a safe starting guess for the Newton-Raphson solution of Eqs. 
(|ITD-(|r8D in the ordered phase, due to the fact that some solutions may suddenly 
disappear. 
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Since the transition from anomalous to normal hysteresis in going from the idr to the ndr 
is preceded by the disappearance of a pair of solutions with S m = 0, the line in the phase 
diagram at which these disappear (dashed line in Fig. [I]) provides an estimation of the place 
at which the former transition occurs. Of course both phenomena are different and so the 
disappearance of a pair of solutions with S m = does not imply an anomalous-to-normal 
transition (recall what happens at the left of the dip in the dashed line of Fig. |l|). 

Admittedly, all of our results are mean-field ones. Although this approximation shows 
undoubtedly its ability to reveal the richness of the phase diagram of this model, it is 
reassuring to see that those of our results that are not original do coincide with the numerical 
simulations for the anomalous hysteresis loop shown in Ref. ||. Nonetheless, the ultimate 
verification of these amazing and potentially useful phenomena lies on the experimentalists' 
side. We hope to see advances in that direction in a near future. 

ACKNOWLEDGMENTS 

The authors thank V. Grunfeld for a critical revision of the manuscript. Partial support 
for this work was provided by the Argentine agencies CONICET (grant PIP 4953/97) and 
ANPCyT (grant 03-00000-00988). 



VI. APPENDIX 



A. The stationary probability distribution function 

We shall adopt as our standard reference the book by Risken |§, whose equation (3.67) 
we have written in the form of Eq. fl7|), with R(x) and S(x) defined in Eqs. At 
variance with Risken's choice (and in accord with Ref. |J) we have adopted q — 1 in our 
Eq. (H), equivalent to its Eq. (3.68). Hence Eqs. (3.95) are translated into 

DW(x) = R+^SS', (23) 
D^(x) = 1 -S 2 , (24) 
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(the prime indicates a derivative with respect to x) whereupon (D^)' = S S'. 
The Fokker-Planck equation (4.46) can be written as 



d t P{x,t) = -d x J{x,t) 



(25) 



with 



J(x, t) = D w (x) P(x, t) - d x [D {2) (x) P{x, t) 



R+^SS'^J P{x, t)-SS' P(x, t)~\s 2 d x P{x, t) 
1 



= RP(x,t)--Sd x [SP(x,t)]. 
On account of Eq. (^), the stationary case of Eq. (|25j ) (namely dtP 



x 



J(x,t) = const. = J, whence 



2J 



s 2 



P st {x). 



(26) 
0) implies 

(27) 



This equation has the form 



y'(x) = a(x) + f3(x)y(x), 



(28) 



and its general solution is 

(29) 

As long as a(x) ^ 0, the integration constant K can be chosen so that y(0) = 0. Otherwise 
(in our case for J = 0) the solution is 













exp 


[ X P(x')dx' 


| J dx"a(x") exp 


- r f3(x')dx' 






Jo 




Jo 





y(x) = K exp 



f3{x')dx' 



If a(x + L) = a(x) and (3{x + L) — f3(x) then 



X+L p(x')dx' = [ L ' f3{x')dx + f X+L P(x')dx', 
o Jo Jl 



x+L 



dx"a(x") exp 



f3(x')dx' 



px rx" 

= / dx"a(x") exp — / f3{x')dx' 
Jo ' Jo 

px-\- L px" 

+ / dx"a(x") exp — / (3(x')d. 
Jx Jo 



(30) 



(31) 



(32) 
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The first term in Eq. fl3T|) is just a number N and the second equals Jq /3(x')dx' (as can be 
seen by writing x' = x" + L). Hence 

i-x+L 



y(x + L) = e N \ y{x) + exp 



/ f3{x')dx' (* dx"a{x")exp - f" (3{x')dx 

JO J Jx JO 



(33) 



By imposing the boundary condition y(x + L) = y(x): 

iN r rx if r x+L 



y(x) 



T^ exp 



f3(x')dx' 



dx"a(x") exp 



(3(x')dx' 



(34) 



For our case it is a(x) = -2J/S 2 , (3(x) = (2R/S 2 ) - (In S)' and y = P st whence, by 
calling 



(x) = -2 / „}, L dx' 



o S 2 (x'^ 



N = -<f>(L) holds, and Eq. fl34j) reads 
P st (x) -- 



2Je-*M ( exp\-d)(x)} r+ L , exp 

dy 



1 _ e -<P(L) 



S(x) 



J is fixed by normalization 
J = 



1 _ e <KL) 



9 T L / 2 r?-r- eX P[~^W1 / C X+L ^, cx P^fa)] 1 

Z J-L/2« X S(a;) \Jx ; 



S(y) 



1 _ e ^( L ) 
2Af ' 



(35) 



(36) 



(37) 



B. Contribution to 0(-L) of the interparticle interaction 

For < A < 1/4, in Eq. (|4]) has a minimum height —(1 + A) at x = and 

a maximum height 1 — A at x = ±7r. The critical points of its derivative are shifted an 
amount of order A towards x = from their A = position of x = ±vr/2, and so S 2 (x) 
has its two maxima (of height 1 + Q/T) inside the cosx > region. Now [S^x)] -2 has two 
minima of height (1 + Q/T)~ l in that region, whereas in the cos a; < one it remains of 
order one. 
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C. Small-x expansion of (j)(x) 



To first order in y, Eqs. (|) and (JTI| read = + K C m + 4A) + K ^m + F, 

S 2 ^) = 2T[1 + (Q/T)(l + 44) V]- Hence in this approximation 

^(*) = r(1 ^^ + ^f m) ln[l + (Q/T)(l + 44) V] 



2Q(l + 4Ay 
(K S m + F) 



T/Q 



arctan[WQ/T(i + 44)] 



1 + 44) 

(1 + 44 + K C m )x 2 /2 - (K S rn + F)x. 



(38) 



Assuming K S m = F = we may approximate 



P st ( x ) = _L 



1 exp[— 0(x)] 

A7 



V / 2TA/^ 



■. exp 



;i + 4A + x c m )x 2 

2T 



1_^_(1 + 4A)V 



(39) 



which clearly has a maximum at x = 0. For small F the maximum shifts toward (K S m + 
F)/(l + 4v4 + KoC m ) and that will in turn produce a small shift in S m from zero, in the 
direction that the maximum shifts (i.e., that of F) whose consequence (by the argument in 
the last subsection) is a reversed current. 

D. Calculation of the particle current 



According to Eq. (3.85) of Ref. @, 



D^(x,t) = hm-(X(t + T)-X(t)) 



X(t)=x 



(*(*)> 



X(t)=x 



(40) 



(in our case, since D^\x) does not depend explicitly on time, this is the conditional average 
over realizations of the noise for any t). The unrestricted average at time t is then 

i-L/2 



(X(t)} - 
From Eq. (p6j ) it is 



-L/2 



dxD {1 \x)P{x,t) 



[L/2 




/ dx 


R + -SS' 


l-L/2 





P(x,t). 



(41) 



(X(t)) = f L/2 dxJ(x,t) + [L> (2) (x)P(x,t) 

J-L/2 L 



L/2 
-L/2 



(42) 
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so for periodic P(x, i) it is 

fL/2 

(X(t)) = / dxJ(x,t) (43) 

J-L/2 

and in the stationary state, where J(x,t) = const. = J: 

(X) = JL. (44) 
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FIGURES 

FIG. 1. Phase diagram of the model for T = 2, A = 0.15 and F = 0. The ordered region lies 
above the full line. Above the dashed line there may exist up to three solutions when S m ^ 0, 
whereas below it there may exist at most one. The squares represent states at which the shape 
of P st {x) and the behavior of (X) as a function of F have been investigated. They correspond to 
Kq = 10 and Q =1, 3, 6, 9, 12, 16 and 21 respectively. 

FIG. 2. Illustration of the passage from the idr to the ndr, as the noise intensity Q increases 
for Kq = 10 (solid line: Q = 5.95; dashed line: Q = 6.5). Only the solution with C m < survives 
after the disappearance of the other two. 

FIG. 3. P st (x) for K = 10 and the values of Q in Fig. 1. For Q = 1 and Q = 21 (solid line) it 
is symmetric, being asymmetric for the remaining values (dashed lines). For Q between 6 and 9 it 
becomes bimodal and as Q increases the peak with larger |x| overtakes the other one, although it 
never reaches beyond 0.5 

FIG. 4. The order parameter V m = (X) (particle current) as a function of Q for Kq = 10 
and F =0 (squares), 0.3 (upward triangles) and 0.44 (downward triangles). The vertical thick line 
signals the transition from anomalous to normal hysteresis. 

FIG. 5. (a) V m = (X) vs F = for K = 10 and Q = 5.97 (just on the left of the dashed line 
of Fig. 1). The stable solutions are those with larger V m values; the other two solutions lie on the 
branches of the C m > and C m < ones for S m = (the solutions with C m ps are not included) . 
(b) Same as for Q = 6.0: one of the unstable solutions has receded from the F k, region (together 
with the C m w one, not shown), (c) Same as for Q = 6.1, showing a complete recession from the 
F pt region, (d) Same as for Q = 6.5: not until the dotted line has completely disappeared do 
solutions in the stable branch appear for |F| > 0.5 and normal hysteresis sets in. 

FIG. 6. V m = (X) vs F for Kq = 10 and Q =1, 1.7, 3, illustrating the appearance of multiple 
solutions and of anomalous hysteresis as the ordered region is reached from the left. 
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FIG. 7. V m = (X) vs F for K = 10 and Q =9, 12, 16, 21, illustrating the disappearance of 
the normal hysteresis loop at the rentrance (the disordered region is reached from the left). 

FIG. 8. (a) V m = (X) vs F = for Q = 6.0 and K =7.25, 8.0 and 10.0, illustrating the way the 
transition from normal to anomalous hysteresis proceeds as Kq increases at Q = const, above the 
dip in Fig. 1. (b) Same as for Kq =14.0 and 18.0: the cusp at the endpoints of the anomalous loop 
has developed into a curl, thus reducing further its excursion. 

FIG. 9. V m = (X) vs F = for Q = 9.5 and K =16.0 (idr) and 13.0 (ndr). For K = 13.0 there 
is still a remnant of the curl existing in the anomalous zone. 

FIG. 10. Illustration of the passage from idr to ndr, as the coupling Kq decreases for Q = 10. 
As in Fig. 2, only the solution with C m < survives after the disappearance of the other two. 

FIG. 11. V m = (X) vs F = for Q = 4.0 (at the left of the dip in Fig. 1) and K =9.0 (idr) and 
8.2 (ndr). One of the branches of unstable solutions has disappeared (together with the C m « 
one, not shown), but the hysteresis loop remains anomalous. 
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